High Jet Efficiency and Simulations of Black Hole Magnetospheres 

Brian Punsly^ 



ABSTRACT 

This article reports on a growing body of observational evidence that many 
powerful lobe dominated (FR II) radio sources likely have jets with high efficiency. 
This study extends the maximum efficiency line (jet power fa 25 times the thermal 
luminosity) defined in Fernandes et (2010) so as to span four decades of jet 
power. The fact that this line extends over the full span of FR II radio power 
is a strong indication that this is a fundamental property of jet production that 
is independent of accretion power. This is a valuable constraint for theorists. 
For example, the currently popular "no net flux" numerical models of black 
hole accretion produce jets that are 2 to 3 orders of magnitude too weak to be 
consistent with sources near maximum efficiency. 

Subject headings: Black hole physics — magnetohydrodynamics (MHD) — galax- 
ies: jets — galaxies: active — accretion, accretion disks 



Introduction 



Relativistic jets emanating from AGN (active galactic nuclei) exist in a variety of 
strengths. Most quasars are radio quiet objects that have either no measurable jets or jets 
that are so weak that they often cannot propagate out of the host galaxy. About > 10% of 
AGN have highly luminous radio jets of which 20% are classic FR II radio sources defined 
by jets that propagate hun dreds of kpc, terminating in lobes of plasma with similar linear 
extent (IdeVries et al.ll2006l ). The energy flux in these jets, Q, can be enormo us with many 



independent estima tes finding long term time averages, Q > 10 ergs/sec (IWillott et al. 



19991 : iPunslyl l2007d ). These jets are not perfectly steady, so there are episodes in which 
the instantaneous power, Q{t), must be even larger. In this Letter, an attempt is made to 
expand on and consolidate the evidence for large jet power that has been accumulating in 
the literature since 2006. 
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The second section quantifies Q relative to the dynamics of tlie accreting gas. From 
observations, one can estimate tfie ratio of Q to the thermal luminosity of the accretion flow, 
Lfoo/, defined as i? = Q/L}yoi or with respect to the mass accretion rate Q = tiqMc^, where 
7] is called the efficiency. Strong radio source s can be described in terms of a concept of 



maximum jet efficiency that was introduced by lFernandes et all (120101 ). i.e a maximum value 
of R or rjQ for radio jets. The maximum jet efficiencies implied by various lines of research 
are compiled. The results are analyzed for consistency and are critically examined in terms 
of the assumptions and limitations of each method. 

In the third section, we discuss the high efficiency sources from section 2 in the context 
of current numerical work, the 3-D numerical simulations of MHD (magnetohydrodynamic) 
accretion flows around black holes with no net magnetic flux. The nexus between these 
simulations and observation is that the simulated Q is always expressible in units of Mc^ 
by all research groups. By making reasonable assumptions about rith, the thermal efficiency 
{Lboi = rjthMc^), as derived by accretion disk theory and new turbulent MHD simulations, 
one can compare the observations with the simulations. 



2. Evidence for High Jet Efficiency 



Ever since the seminal work of iRawlings and Sanders I ( ll99ll ). astrophysicists have been 
trying to estimate the enormous energy flux that feeds the radio lobes in FR II radio sources 
and relate it to the thermal luminosity of the accretion flow. The three most viable options 
for estimating Q are eit her based on the low frequency (151 MHz) flux from the radio lobes 
on 100 kpc scales (eg. IWillott et al.l (Il999l )). or the work done creating the cavities that 
are ca r ved out of the in t ra-clu ster medium by the expanding radio lobes (eg. iBirzan et al 
( 120041 ): iMcNamara et all (|2010[ )). or models of the broadband Doppler boosted synchrotron 
and inverse Co r apton radiation spectra associated with the relativistic parsec scale jet (eg. 
Ghisellini et all ( 120101 )). Each method has its advantage. The 151 MHz method is the most 
widely applicable, all that is needed is a radio spectrum. A disadvantage is that it involves 
long term time averages, Q, that do not necessarily reflect the current state of quasar activity. 
The second method is also not contemporaneous, yet more accurate in principle than the 
flrst, but is restricted to low redshift sources with deep X-ray observations. The last method 
is contemporaneous, so one can deflne a directly interpretable ratio of jet to accretion thermal 
power, R{t) = Q{t)/Lboi- However, one is forced to deal with estimating the large Doppler 
enhancement factor, 6, which is a potential source of large uncertainty (the jet luminosity 
scales like 5^). In this section, usi ng these three method s, we expand on the notion of a 
maximum jet efficiency deflned in iFernandes et all (120101 ) that is based on 151 MHz flux 
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estimate. 



In Figure 1, the black squares are a scatter plot, Q versus L h „i, of the complete sample 
of FRII narrow line radio galaxies (NLRGs) from iFernandes et all (|2010[ ). The low frequency 
selected sample is limited to0.9<2<l.l,asa compromise between having sufficient cosmic 
volume to find strong radio sources and being sufficiently close so that these sources can be 
detected in the IR. L m. = 8.5z/Lj,(12/im) is computed from the IR luminosity at 12 /im 
Richards et al I (120061 ). They define a diagonal line at i? = Q/Lhoi = 25 (the black solid 



line in Figure 1), the maximum efficiency line. The dashed blue vertical line represents the 
approximate dividing line between Seyfert 1 galaxy and quasar luminosity [My = —23). The 
dashed vertical orange line represents the dividing line between solely Seyfert 1 galaxies and 
a mixture of LLAGNs (low luminosity AGN) and weaker Seyfert 1 galaxies (iHol l2005l ). In 
contrast to Seyfert 1 galaxies and quasars, the LLAGNs do not have a strong "bl ue bump" 



in their spectra, the si gnature of strong thermal emission from an accretion flow (IHol 12005 



Sun and Malkanlll989l ). The LLAGNs are estimated to have inefficient modes of accretion, 



expressed in terms of the Eddington luminosity as Lboi/LEdd < 10 ^ in contrast to quasars 



and Seyfert galaxies which typically have 10 ^ < Lhr^i/L pM < 1 (|Holl2005t ISun and Malkan 



19891 ). Thus, Figure 1 indicates that the NLRGs in the IFernandes et all (120101 ) sample are 
likely to have a central black hole in a high efficiency accretion state typical of a quasar or 
Seyfert 1 galaxy, but the optical/UV core is hidden. 



One can expand the IFernandes et all ( 120101 ) treatment to a larger range of black hole 
accretion states and jet power, b y considering the low redshift sample of IR observations of 
FRII NLRGs oflOgle et all (120061). T he same IR bolometric correction and Q estimators can 
be used as in IFernandes et all (120101 ). They noticed a diagonal boundary in a scatter plot 
of IR luminosity versus 178 MHz flux which is similar to the maximum efficiency li ne. I n 
Figure 1, the orange circles represent these "weak Mid-IR sources" from I Ogle et all ( 120061 ). 
Notice how well the orange circles respect the maximum efficiency line. The trend now 
extends below the upper limit of LLAGN luminosity. Since Seyfert 1 galaxies also exist at 
such luminosity, and the trend looks smooth, there does not seem to be any evidence of a 
change in accretion mode for FRII NLRGs at Lboi < 6 x 10^^ ergs/sec. 

A small sample of FRII NLRGs from McNamara et all ( 2010 ) is also plotted as blue 
triangles in Figure 1. In this sample, Q is estimated by an independent method, the work 
done to create la rge bubbles in t he intra-cluster medium. The IR luminosity is estimated 
from the data in IShi et al I (l2005l) with the sync h rotro n component subtracted off and the 



IR bolometric correction is from iRichards et al 
maximum efficiency line concept. 



([20061). This data also conforms with the 



Another method of quantifying a jet as highly efficient is to choose sources with QEdd = 
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Fig. 1. — The maximum efficiency line is illustrated in this scatter plot of Q and Li,oi- The 
scatter plot shows that 1 < i? < 25 is not an unusual state of jet activity. 
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Q/LEdd > 1- This is an extreme conditi on since Lm/Lem > 1 quasars ar e either extremely 



rare or as is often argued nonexistent ( iMarconi et all l2009l : iNetzed 120091 ). Thus, the Qsdd 
sources would almost certainl y have instanta neous episodes with R{t) > 1. A small sample 
of these sources were foun d in iPunslv J2007ch. These are the red diamonds in Figure 1 and 
overlap the high end of the iFernandes et all (120101 ) scatter. They are lobe dominated quasars 
for which UV continuum emission and broad line st rengths were used to e stimate Lh„i an d 



151 MH z flux was used to estimate Q. Note that the IFernandes et all (l2010l ) and the IPunsly 
(l2007cl ) samples are consistent with a maximum value of Q > 10^'' ergs/sec that seems to 
make the maximum efficiency line bend over towards the equipartition line defined by i? = 1. 

The fact that Q seems to reach a maximum does not necessarily mean that there are not 
even larger instantan eous iet powers . By fi tting broadband blazar spectra, from the radio 
band to gamma rays, iGhisellini et all (120101 ) believe that they have a method to extract the 
jet power within a few light years of the central black hole in a blazar - almost contempora- 
neous on cosmic time scales. The model is one of a highly relativistic magnetized plasmoid 
propagating in the radiation environment of the quasar. The inverse Compton emission is 
necessarily modeled simultaneously with an accretion disk model of the "big blue bump" for 
each source. The method is a bit controversial because of the large Doppler enhancement in 
blazar jets and the uncertainty that it introduces in the intrinsic luminosity. The appeal of 
this method is that it is completely independent of the techniques used for the other data 
in Figure 1. In order not to clutter Figure 1, only the R{t) > 2 blazars are plotted. These 
Q{t) estimates (that are plotted as dark blue circles) also respect the maximum efficiency 
line with only two outliers. 

None of the methods used to create the data sets in Figure 1 is a rigorous justification 
of the maximum efficiency line in isolation. However, the agreement that is achieved by 
these independent experiments are strong scientific eviden ce in support of the notion of the 
maximum efficiency line found in IFernandes et all (l2010l ) that is now extended to over 4 
decades in Q. The fact that this line extends over the full span of FR II radio power is 
an indication that this is a fundamental property of jet production that is independent of 
accretion power. Another important aspect of Figure 1 is that powerful FR II radio sources 
are plentiful in the range 1 < i? < 25. Furthermore, since jet power is not steady over the 
lifetime of the QSO, many of the sources below the R = 1 line likely have episodes in the 
high efficiency range 1 < i? < 25, i.e. this is not an aberrant or outlier state of jet activity. 
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~ Doubly Truncated Ergospheric Disk 
° KDJ Ergospheric Disk 

McKinney Blandford-Znajek 

' Hawley et al Blandford - Znajek 
O Hawley et al Counter-Rotating Disk 



• KDE Ergospheric Disk 
' KDH Ergospheric Disk 
" ~ Equipartion Line 

" " "Ergospheric Disk: Inner Boundary Normalized 
Maximum Jet Efficiency 



Fig. 2. — A comparison of 3-D numerically simulated data with the constraints on jet 
power in terms of accretion rate (Q/Mc^) imposed by observations. The "region of high jet 
efficiency" corresponds to the region between R = 1 and R = 25 in Figure 1. 
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Theoretical Discussion 



In the past, we were forced to compare a scatter plot like Figure 1 to theory by means 
of parametric niodels. However, these models have a l a rge unknown, the st rength of the 



magnetic field (IBlandford and Znajekl 119771 : iMeierl l200ll : iNemmen et all |2007| ) . The ability 



versial, rendering the flux distribution as a major unknown ( 


Ghosh and Abramowicz 


1997; 


Livio et al 


1999 


Rothstein and Lovelace 


2008; 


Revnolds et al 


2006). Long term MHD sim- 



ulations in a generally relativistic background can at least provide a self-consistent magnetic 
field distribution (be it not necessarily unique). Thus, the best tools that we have to inves- 
tigate the central engine of radio loud AGN, without the over-riding uncertainty of the field 
distribution, are the current battery of long term 3-D numerical simulat i ons of black hole ac- 



2006; Krohk et al 



2005 



cretion systems (McKinnev and Blandford II2009I ; iBeckwith et alll2008bl ; iHawley and Krolik 



Fragile et all 120071 ). The initial state is a thick torus of gas in equi- 
librium that is threaded by concentric loops of weak magnetic flux that foliate the surfaces 
of constant pressure. There is no net magnetic flux in these simulations. If the loops are 
configured in the same orientation and are poloidal and not toroidal then t he leading edge 
that a ccretes will deposit a net poloidal flux on the black hole and a jet forms ( iBeckwith et al 



2008bl ). In these simulations, Q is expressible in terms of accretion rate onto the black hole. 



Mc^. Thus, one can compare different simulations a nd one can compare to the o bservations 



if rjth is known. Theoretically, r]th was determined by lNovikov and Thornel ( 1l973l ). This was 
recently ques tioned since the boundary condit ion of zero stress at the inner most stable orbit 
was suspect ( lGammielll999l ; iKrolik et allll999f). In spite of this, recent high resolution simu- 



lations of accretion disks indicate that the iNovikov and Thornel ( 119731 ) value is accurate to 
within a few percent for thin disks and only increases modestly for high spin rates (expressed 
as a/M < 1, where "M" an d " a" are the black hole mass and specifi c angular momentum, re 



fectively) and thick disks ( IBeckwith et alll2008al : iPenna et alll2010l; iNob 
Therefore, in the context of this discussion, the iNovikov and Thorne 



e and Krolikll2009r) . 



(119731 ) value of i]th 



is suitable for the purpose of comparing the simulations to the maximum efficiency line in 
Figure 2. Recall that the sources in Figure 1 have Li^oi consistent with large viscous dissipa- 
tion (big blue bump), so a thermally lum inous disk model is app ropriate as opposed to low 
efficiency advection dominated accretion (iNarayan and Yi Ill995l ). 



Figure 2 compares the 3-D simulated data to the maximum efficie ncy line. ThelMcKinney and Blandfoi 



( 120091 ) a/M=0.92 simulation shows a Blandford-Znajek (B-Z) jet as in IBlandford and Znajek 



^Using the larger simulated rjth values at high spin rates will just elevate the equipartition and maximum 
jet efficiency lines in the plot slightly, even farther from the already nonconforming simulation data. 
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(1977). They note that Q ^ O.OlMc^ is similar to the 2-D solutions re ported in iMcKinnev 
(l2005l ). Thus, the spin dependent efficiency equation, equation (3) of iMcKinneyl (120051 ) is 
plotted above. The other B-Z je t data (collectively r e ferred t o as Hawley et al. Blandford- 
Znajek in Figure 2) comes from iHawley and Krolik I (120061 ): iKrolik et all (120051 ) except for 
the a/M = 0.99 and the a/M=0.998 dat a. The reason for t he new d ata p oints is that the 
high spin ca ses a/M=0.998 (KD E) form iKrolik et all J2005h: IPunslvl (l2006h and a/M=0.99 
(KDJ) from Iflawlev and Krolikl kood \: IPunslvl J2007ah : IPunslv et a? toidj have a strori g 
ergospheric disk jet (as indicated in the Figure2) that suppresses the B-Z jet (|Punslyll2007bl ). 
To find the p ower of a pure B-Z jet, the raw data from two simulations (that were used for 



ray tracing in lBeckwith et all (l2008al )) that was generously provided to this author by John 



Hawley, KDEb (a/M=0.998) and KDJd (a/M=0.99), was reduced. These simulations were 
different from KDE and KDJ because the code was modified to include artificial diffusion 
terms in th e equations of con tinuity, energy conservation, and momentum conservation as 
described in lDe Villierd (120061 ). The resultant numerical diffusion suppressed the ergospheric 
disk jet giving a more pristine estimate of the B-Z efficiency than can be obtained from KDE 
and KDJ. The simulation, KDH (a/M=0.9 5), has a weak ergospheric disk jet so the esti- 
mates for the B-Z power are straightforward iPunslyl (l2007bl ). The mechanical energy flux is 
i ncluded in the se estimates which can be non- negligible in some of the high spin simulations 
(IPunslyl l2007bl ) . Note that no accretion disk jets form in any of these simulations. 



Notice that all the no net flux simulations in Figure 2 fall below the region of high jet 
efficiency, 1 < i? < 25 from Figure 1, being 2 to 3 orders of magnitude less efficient than the 
most efficient jets. The only simulations that are close are the high spin ergospheric disk 
jets from KDE and KDJ. But they are still just below the region of high jet efficiency. Can 
the a/M=0.998 case be optimized to reproduce the high jet efficiency? To answer this, we 
explore the curious situation that the ergo spheric disk jet in KDE is not much stronger than 
KDJ contrary to the theory (IPunslyl l2008l ). A spin dependent expression for the ergospheric 
disk jet Poynting flux (approximately Q in this discussion), S, can be approximated as 

S ^ N{B') [[SA]Qh]' , (1) 

where the vertical magnetic fleld in the ergospheric disk is B, the surface area of the ergo- 
spheric disk is SA and the fleld line angular vel ocity scales w ith the horizon angular velocity, 
Qh, and N{B'^) is a normalization constant (IPunslyl l2008l ). The theory does not flx the 
distribution on B in time and space and SA of the jet producing region that will occur in 
a given numerical simulation. In equation (1), the empirical distribution of B^ is absorbed 
in the normalization constant A^. In principle, SA of the equatorial plane in the ergosphere, 
changes very rapidly with spin f or a/M < 1, hence the expectation of signiflcantly larger jet 
power for KDE than for KDJ ( IPunslyl l2008l ). Setting N = 0.0028 in equation (1) creates 
the "doubly truncated ergospheric disk," a two parameter {SA and Qh) fit to the simulated 
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data in Figure 2. It is doubly truncated because the jet does not fill the entire ergospheric 
equatorial plane, but is restricted to the region rj„ < r < 1.50M, where rj„ is the inner calcu- 
lational boundary (located outside the event horizon, r^n > th). For some unknown reason, 
the vertical flux that creates the jet dies off rapidly beyond r = 1.50M in the simulations. 
Yet, in princi ple, an ergospheric disk can exist throughout the ergosphere th < r < 2M 
( lPunslyll2008l ). The simple two parameter equation (1) (with N fixed) is a reasonable fit to 
the numerical data (the 3 blue dashes in Figure 2) and explains the dependence of S on a. 
Equation (1) implies that the reason KDE is not much stronger than KDJ is because SA in 
the computational grids of the two simulations is virtually identical (~ 18M^ in geometrized 
units). If one were to move rj„ in KDE inward so that it is normalized to the the KDJ 
scaling of rin{KDJ) = 1.203M = 1.054r//, one expects a larger SA and therefore a larger 
jet power, (i.e. change the inner calculational boundary in the numerical grid of KDE from 
rin{KDE) = 1.175M to r^^ = l.ObArn = 1.12M). SA is increased to 23.3M^ for a/M=0.998 
if Tin = 1.12M < r < 1.5M. The black dashed curve in Figure 2 is a plot of the predicted 5* 
from equation (1) with the empirically fit value of = 0.0028 from the simulated data and 
normalized numerical grids defined by rj„ = l.ObArn- Even with the enhancement in power 
for a/M = 0.998, with a renormalized grid, the curve barely penetrates the region of high jet 
efficiency. Based on this analysis of this most efficiently known simulated jet, it is concluded 
that any attempt to optimize the no net flux scenario will still be too weak to explain the 
region of high jet efficiency in Figure 2. 



4. Conclusion 



This study shows that the maximum jet efficiency condition, R ^ 25, extends over the 
entire range of known FR II jet powers. It is demonstrated that the no net flux accretion 
numerical models can not explain the large number of high efficiency jets with R > 1. 
Furthermore, the initial no net flux state in the considered simulations is configured to 
yield the maximum jet power. The situation for no net flux is even more nonconforming 
if the initial field is not composed of loops of the same vertical orientation, but randomly 
oriented. In this case, the jet power can be reduced by three orders of magnitude or more 
( iBeckwith et alll2008bl ). It is concluded that the no net flux models are not suitable numerical 
models for powerful radio loud quasars. Clearly some other dynamical element is needed in 
the setup of the sir nulations. Perh a ps it is the accretion of a large reservoir of large scale 
magnetic flux as in llgumenshchevi (120081 ). This ca n have two releva nt effects. First, the 
ergospheric disk power should be greatly increased (IPunsly et all 120101). Also, u nlike the no 
net flux simulations, a magnetized accretion disk forms in llgumenshchevi ( 120081 ). 
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Finally, note that Figure 2 does not indicate that large "a" acts likes a switch that 
transforms a radio quiet black hole accretion system into one with high jet efficiency. The 
only thing that looks remotely like a switch in Fi gure 2 is the ergospheric d isk. This con- 



clusion is in accord with the B-Z switch model of iTchekhovskoy et all (|2010[ ) that requires 
radio quiet (loud) quasars to have spins a ~ 0.15 (a ~ 1). This scenario seems unlikely, ra- 
dio quiet quasars are typically selected by optical/UV luminosity and therefore shoul d have 



elevated (mass and angular momenturn) accr etion rates and large a (IBardeen Ill970l ). The 



magnetospheres in ITchekhovskoy et all (120101 ) are shaped by a boundary condition. There 



is no accretion in the simulation, therefore there is no calibrated measure of the jet strength 
with respect to Mc^. 

I would like to thank an anonymous referee who offered many constructive comments. 
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